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Abstract 

l-H , 

C/3 

CJ I The compute-and-forward (CMF) method has shovi^n a great promise as an 

^\j ' innovative approach to exploit interference to^vard achieving higher net^vork 

OO ■ throughput. The CMF was primarily introduced by means of information 

<N . 

^T . theory tools. While there have been some recent works discussing differ- 

^. ' 

Tlj" , ent aspects of efficient and practical implementation of CMF, there are still 

o 

f^ I some issues that are not covered. In this paper, w^e first introduce a method to 



decrease the implementation complexity of the CMF method. We then evalu- 
ate the exact outage probability of our proposed simplified CMF scheme, and 
hereby provide an upper bound on the outage probability of the optimum 
CMF in all SNR values, and a close approximation of its outage probability 
in low SNR regimes. We also evaluate the effect of the channel estimation 
error (CEE) on the performance of both optimum and our proposed simpli- 
fied CMF by simulations. Our simulation results indicate that the proposed 
method is more robust against CEE than the optimum CMF method for the 
examples considered. 

Index Term- compute-and-forward, outage probability, performance anal- 
ysis, channel estimation error, wireless relay network. 
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I INTRODUCTION 

The CMF method, proposed by Nazer and Gastpar [1], enables exploiting, rather than 
combating, the multiple access interference in a wireless relay network, and thus results in 
improved network throughput [2]. In this method, relays, instead of recovering single messages, 
attempt to reliably recover and pass an integer linear combination of transmitted messages, 
called an equation, to the destination. By receiving enough equations, the destination can 
solve the hnear equation system to recover desired messages. 

In recent years, the CMF method as a promising approach, has received a lot of attention. 
There have been an increasing number of works on theoretical aspects of CMF method. In [3], 
by analyzing asymptotic behavior, the number of degrees of freedom of CMF method is derived. 
The authors in [4] generalize the CMF method to the case of multiple antenna sources and 
relays, and compute the corresponding achievable rate based on an optimization problem. 
Successive recovering of messages in relays for CMF method is proposed in [5]. The idea is 
similar to successive interference cancellation (SIC) technique and leads to a higher recovering 
rate at each relay. In [6], the impact of inter-symbol interference on the CMF method, for the 
bi-directional relay case, is evaluated. 

Most of the analytical studies, including aforementioned works, consider the CMF from 
an information theory viewpoint which involves finding an achievable rate or a capacity re- 
gion [7-9]. On the other hand, some works concentrate on the practical aspects and implemen- 
tation of the CMF [10]. Efficient design of proper codes and lattices are discussed in [11-14]. 
Reference [15] considers the problem of lattice decoding and proposes some practical and effi- 
cient approaches to this end. Finding an optimum integer Equation Coefficient Vector (ECV), 
coefficients of recovered combination at a relay, is computationally complex. This issue in not 
addressed considerably and will be investigated in this paper. 

The authors in [1] and [16] have calculated the outage performance of the CMF method in 
the cases of three-transmitter multiple-access channel (MAC) and Multi-Way relay channels, 
respectively. However, in these works, the achievable rate of CMF has not been given by a 
closed-form expression, but rather as the solution of an optimization problem [1]. In fact, 
to the best of our knowledge, all the previous works, including [1] and [16], have employed 
numerical calculations or simulations as a part of their analysis to find the outage probability. 



In this paper, we first present some lemmas which help us to reduce the complexity of 
optimal CMF method. These lemmas are then used to approximate the optimum CMF as a 
simplified CMF with considerably less computational complexity. We provide the exact outage 
probability of our proposed scheme and show that our proposed method performs near the 
optimum CMF in low SNR values. We also consider the effect of the channel estimation error 
(CEE) on the performance of both the optimum and simplified CMF methods, our simulation 
results show that the proposed scheme is more robust against the CEE than the optimum 
scheme for the examples considered. 

The rest of this paper is organized as follows. Section 11 introduces the system model and 
assumptions. The optimum and the proposed CMF methods are presented in section 111. Sec- 
tion IV includes the performance analysis of the proposed method. Simulation and numerical 
results are presented in section V. Finally, the paper is concluded in section VI. 

II SYSTEM MODEL 

We consider a general network, which itself can be a part of a larger network, consisting of 
L = 2 transmitters as message sources, M relays, and one receiver (destination), as shown in 
Fig. 1. The destination aims to reliably recover both messages of the sources. 

Let the channel gain from source / to relay m and the channel gain vector of relay m be 
denoted by hmi and h^ = [hmi,hm2]'^, respectively. We consider block fading channels, and 
assume that the channel gains are real independent and identically Rayleigh distributed (i.i.d.) 
variables with unit variance. The channel noises are additive white Gaussian (AWGN) with 
unit variance. Let Pi be the transmission power of source /. Then, the average SNR received at 
each relay from source / is equal to Pi, when not considering the shadowing effects or assuming 
the same shadowing effects for all relays. The results can be easily extended to more general 
cases. The channels from relays to the destination are assumed to be orthogonal and have high 
enough capacity (i?o in Fig. 1) to reliably transfer the required information. 

Our proposed method of transmission over the network is based on the CMF scheme. A 
detailed description of CMF scheme can be found in [1]. First, two sources map their messages 
wi and W2 to symbols xi and X2, respectively, and transmit the symbols simultaneously. 
Receiving a noisy linear combination of transmitted symbols y™,, the m-th relay, m = 1,2, , M, 



computes an equation (integer linear combination) of transmitted messages with ECV equal 
to S-m = [01,02]^ G Z^.Then, all relays pass their computed equations and ECVs to the 
destination. Finally, the destination attempts to extract both messages from M (M > 2) 
received equations. Specifically, the destination selects two independent equations with the 
highest rates and solves them to recover both the messages. 

Ill EFFICIENT ECV 

We attempt to find an efficient ECV in the relay that provides the highest computation 
rate (the rate of recovering an equation) while imposes tolerable complexity on the system. A 
relay m with channel vector h^ can recover an equation with ECV equal to a, as long as the 
message rates are less than the computation rate [1] defined as 



R{K 




(1) 



where the two sources have the same power P, i.e. Pi = P2 = P, and 



log^(x) = max(x,0). 



(2) 



If we define the vector g^ as 



[9ml, 9; 
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computation rate for unequal source powers can be written as 



^(gm,a) = -lo, 
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It is notable that g^i is the instantaneous received SNR from source I at the relay and hence, 
||gm|| (called instantaneous sum SNR) equals the sum of all instantaneous received SNRs. In 
the following, the optimum method along with our proposed suboptimum method of selecting 
ECVs is presented. 



A. Optimum Compute-and-Forward 

In the optimum CMF nietliod [1], eacli relay selects the ECV with the highest computation 
rate; More specifically relay m calculates its ECV according to the following maximization 
problem 

Sim = argmaxi?(gm,a), (5) 



where is equivalent to [10] 



a^ = argmina-^Gma, (6) 

agZ2 
a^O 



where gm is defined in (3) and G^ is a positive-definite matrix and is defined as 

GA / T SmSm \ /-x 

V 1+ llgmll / 

The above integer optimization problem is equivalent to the shortest vector problem (SVP), 
and has no closed-form solution [10, 16]. Lattice reduction algorithms such as [17] can be 
applied to calculate the optimum ECV numerically. However, due to the large search space of 
the problem on hand, the complexity of the system is considerably high. In the following, we 
introduce some techniques to reduce the size of search space and then propose our simplified 
compute-and-forward method. 

B. Reducing the search space 

The optimization in (6) is defined over Z^ and has a infinite search space. To limit the 
size of search set, which leads to reduced complexity, we use the following lemmas. Note that 
all the following lemmas are also hold for the more general case when the number of sources 
are larger than 2 (L > 2). In fact for this general case, we have gm = [gmi, ■ ■ ■ , gmi]'^ and 
^m = [(^1, ■ ■ ■ jO-l]^ £ 2'^- Moreover, although we consider real vectors, the complex case can 
be modeled as a real one with 2L sources and 2M relays [1]. 

Lemma 1 (Lemma 1 in [1]) An ECV a^ results in a zero computation rate for 

II 11^ "^^ "1 I 1 1 1 1 -^ / (n\ 

\\^-m\\ ^1+ llgmll • (8) 



Prom this lemma, for the optimization problem of (6), as stated in [1], it is sufficient to 
search only over ECVs that satisfy ||am|| < 1 + ||gm|| • Therefore, using Lemma 1 results in 
a finite search space. 

Lemma 2 For the solution of the optimization problem defined in (6), we have either 

sgn {ami) = sgn (gmi) , for alll = 1,- ■ ■ , L, (9) 

or 

sgn (ami) = - sgn (g^i) , for all I = 1,- ■ ■ , L, (10) 

and if Qmi = 0, then ami = 0. 

The operator sgn(-) shows the sign function. 

Proof The only term in (4) that depends on the sign of elements of a is the term la^gml = 
\Y2i^i9mi\ ■ For maximizing (4), or equivalently (6), when the norm of a is fixed, this term 
has to be maximized. Thus, all the terms aigmi, I = ^r ' ' , L must have the same sign and this 
proves the lemma. | 

Lemma 2 suggests that the search can be done over non-negative integers. Specifically, we 
replace the elements of gm with their absolute values and solve the optimization problem (6) 
by searching over non-negative integers. The solution gives the absolute values of the elements 
of optimum ECV a^, where the signs of its elements are determined according to Lemma 2. It 
is noteworthy that both a^ and — a.^ result in the same computation rates. In the rest of this 
paper, without loss of generality, we assume that all elements of a^ and gm are non-negative. 

Lemma 3 Any ECV a^ with gcd {omi, ■ ■ ■ , cimL) > 1 cannot he the solution of (6). 

The operator gcd (ami, • • • , o-mL) shows the greatest common divisor of integers ami, ■ ■ ■ , cimL 
and we assume gcd(a,0) = a. 

Proof Assume that gcd (ami, ■ ■ ■ , cimi) = g > 1. If we define the integer vector bm = -am, 
then we have 



g2 



and thus a^ cannot be the solution of (6). | 

This lemma limits the search space of (6) as well. It is sufficient to search over ECVs with 

gcd{aml,--- ,araL) = 1- 

Lemma 4 Any ECV Am with ||gmm {^m)\\ > ||gm|| cannot be the solution of (6), where 



rtinm y^m) 



min ||g|| , (12) 

ai;Ga,„<brGb 



l|b|P=£l+||gf 



and 



G^ I--f^ . (13) 

1 + g 



Proof Assume that ECV a^ is the solution of (6). Hence, it satisfies di^Gm^m < b-'^Gmb for 
all b G Z^. From Lemma 1, it is sufficient for a^ to satisfy a^Gmam < b'^Gmb only for all 
||b|| ^ 1 + ||gm|| • Therefore, gm is a feasible point of (12) (since it satisfies both constraints 
of (12)). As a result, for the optimal solution of (12) we have ||gmm (^m)ll < l|gm||- | 

Since ||gm|| is the instantaneous sum SNR. It follows from Lemma 4 that the minimum 
instantaneous sum SNR for which an ECV e^ can be selected as the solution of (6) equals 
||gmm(eA;)|| . Optimization defined in (12) is a mixed integer optimization problem and can 
be solved through numerical methods. Lemma 4 can significantly reduce the search space. To 
show its effect, we have computed the values of ||gmm(eA:)|| for different ECVs e^ up to the 
instantaneous sum SNR about to 2000 (33 dB), in the case of L = 2 sources. The results are 
sorted and provided in Tabel 1. As an example, in Fig. 3, the distance between point F and the 

rp 

origin is equal to ||gmm (efe)|| for e^ = [2, 1] . Now, according to the Lemma 4, for a given g^, 
equivalent to a constant instantaneous sum SNR, only ECVs that satisfy ||gmm (^m)ll ^ l|gm|| 
are needed to be searched. These ECVs can be found from Tabel 1. Note that Tabel 1 is 
computed just once and thus adds no extra complexity to the method. This approach can be 
extended to fading channels, in which the instantaneous sum SNR has a certain probability 
distribution function, and with a high probability it is less than a threshold. This idea is the 
base of our proposed method which is described in the following subsection. 

All of the provided lemmas reduce the complexity of optimum compute-and-forward method 
by limiting the search space of (6). As an example, consider two vectors gi and g2 with 



II gill = 100 (i.e. 20 dB instantaneous sum SNR) and ||g2|| = 1000 (i.e. 30 dB instantaneous 
sum SNR), respectively. Lemma 1 limits the size of search space to 317 and 3141 vectors, 
respectively. Lemma 2 (in conjunction with Lemma 1) decreases it to 89 and 818, respectively. 
The size reduces to 49 and 479, respectively, by using Lemma 3. Finally, by exploiting all of 
the four lemmas, the search set only contains 7 and 23 vectors, respectively. However, as it is 
observed, by increasing the SNR, the search set and complexity grow considerably. 
The following lemma facilitates the completion of Tabel 1. 

Lemma 5 For Any permutation e^ of ECVsk, we have ||gmm(efc)|| = ||gmm(efc)||. 

Proof Let e^ be a permutation of e^, then there exists a permutation matrix U such that e^ = 
Uefe. According to the properties of a permutation matrix, we have that ||efc|| = ||Uefe|| = ||efc|| 
and U^^ = U^. From (12) and (13), we can write 

||gmm(efc)|| = min ||g|| , (14) 

e^Geft<b^Gb 

l|-||2 ,, ,,9 

||b|| ^l+||g||^ 

where 

G^fl--^). (15) 

V 1 + ||g|| / 

By defining g = U^^g = U"^g G M^ and b = U^-'^b = U"^b G Z^, which leads to g = Ug, 

b = Ub, ||g|| = ||g|| and ||b|| = ||b||, we have 



efcGefc = ||efc|| - ^ 

l + l|g|| 

, |g^U^Ue,|^ 

l+l|Ug|| (16) 

17^ 1 2 
-II ||2 |g^fc| 

II fell -, II ii2 

1 + llsll 

= efc'^Gefc 



and, in a similar way, b-^^Gb = b-^Gb. This yields that minimization in (14) is the same as 
the one in (12) and this proves the lemma. | 

As an example, 625 = [2,7] is a permutation of 624 = [7,2] , in Table 1, and hence we 

have llgmm (625)11 = llgmm (624)11- 



C. Simplified Compute-and-Forward (Proposed Scheme) 

As mentioned in the previous subsection, the idea of our method is based on limiting the 
search set of the optimization problem. To describe our simplified scheme, in Fig. 2 we have 
plotted the probability of selecting different ECVs in a relay, employing optimum CMF method, 
versus average SNR , as an example, for the case of equal transmission powers, i.e. Pi = P2 = P, 
in Rayleigh fading channel. From this figure, for SNRs less than 8 dB, the ECV of the optimum 

_^ I rp rp rp ^ 

CMF is in the set < [1, 0] , [0, 1] , [1, 1] > with a high probability. Likewise, for SNRs less than 
16 dB, the ECV of the optimum CMF scheme is in the set |[1, 0]^, [0, 1]^, [1, 1]^, [2, 1]^, [1, 2]^| 
with a high probability. These results are consistent with the results provided in Table 1 as 
well. When the average SNR (of each user) is less than 8 dB, the average instantaneous sum 
SNR, i.e. ||gm|| , is less than 12.7. Thus, according to the Table 1, the first three ECVs can 
be selected as an optimum solution with a high probability. In a similar manner, for average 
SNR (of each user) less than 16 dB, the average instantaneous sum SNR is less than 79.6. 
Thus, based on the Table 1, the first five ECVs cab be selected with a high probability. That 
is Table 1 confirms Fig. 2. 

From the above observations, which is true for the other examples considered as well, we 
propose a simplified CMF, in which the search set for finding the proper values of ECV is finite 
and its size can be properly selected based on the channel received SNR (using Table 1). That 
is 

Arn = arg min a^G^a, (17) 

aGSif 
37^0 

where Sk is the ECV search set including K vectors. For instance, for the example considered, 
we can select 

52 = {ei,e2} , 

53 = {ei,e2,e3}, (18) 
S5 = {ei,e2,e3,e4,e5}, 

and so on, where ei = [1, 0] , e2 = [0, 1] , e3 = [1, 1] , e4 = [2, 1] , and es = [1, 2] are derived 
from Table 1. 

We denote the simplified CMF, with the ECV search set Sk, by CMF{K). In fact, CMF(cx)) 



is equivalent to the optimum CMF. In CMF(i^), the complexity of finding ECV is in the order 
of KM, which is significantly less than that of the optimum CMF. 

It is noteworthy to mention that since CMF(X) limits the search region of ECV selection, 
the outage Probability of CMF[K) upper bounds that of the optimum CMF, for all values of 
K. 

IV PERFORMANCE ANALYSIS 

In this section, we present the exact end-to-end performance analysis of CMF(i^) in terms 
of the outage probability, which provides an upper bound for the outage probabihty of the 
optimum CMF method. 

A. Probability of selecting an ECV in a relay 

From (3) and (6), CMY{K) selects ECV efc, for A; = 1,2, , iC , if and only if i?(g^,efc) ^ 
i?(gm, Gi),'^l 7^ k. Thus the probability of selecting e^ is equal to the Probability that gm is in 
the region D^ defined as 



Dk ^ {g\Rig,ek) ^ Rig,ei),yi ^ k} 

= {g|efc^Gefc < ei^Geiyi ^ k] ,k = I, . . . ,K, 

where G is defined in (13). This probability can be computed as 

Pf ^ Pr {a^ = e4 = Pr {g„ G Du] 

= / f{9mi,gm2) dgmidgm2 ,k = 1,...,K, 



(19) 



(20) 



where f{gmi,gm2) is the joint probability distribution of g^i and gm2- Sine gmi and gm2 are 
independent, we have 

f{gmugm2) = ^-^e-'^.^-^e-'^. (21) 

Region Dk is calculated from (19) by computing the borders e^'^Gek = ei'^Gei.yi ^ k, 
which are in general, in terms of gi and g2, hyperbolas. As an example. Fig. 3 shows the 
selection regions D^ for CMF (5). Note that determining the regions D^, k = 1,- ■ ■ ,K is not 



10 



required for finding a^, just required for the performance analysis. ECV , i.e. a^, is found 
simply by calculating and comparing the terms \^eiJ^GmGk, k = 1,- ■ ■ , K^ according to (17). 

B. The outage probability of each relay 

The outage Probability of a relay is the probability that the relay computation rate, i.e. 
R{gm, ^m), is less than a given target rate Rt, and can be expressed as 



P,% = Pi {R{gm,^m)< Rt}. (22) 



Using the law of total probability, we can write 

K 
P?eTay = ^ ^^^ ^S™ ^ i?fe} . Pr {i?(g„„ a,„) < i?, |g„, £ Dk] 

k=l 



(23) 



fc=i 



where P^^'' is defined in (20) and (a) holds since if gm G Dk then we have a^ = e^. We define 
conditional outage probability as 

Pout\., = Pr {i?(g^, e,) < Rt\g^ e Dk} . (24) 

Using (23) and (24) yields 

K 
rtOut _ \ ^ r>Sel p /oe.\ 



relay / j 

k=l 



By defining region 



Ofc = {g|i?(g,efc)<i?,} 

g|e/Gefe > ^^ } ,k=l,...,K, 



the conditional outage probability Pout\ek can be easily calculated as 

Pout\ek = Pr{gm e Ok\gm ^ Dk} 

Pr {g„ eOkH Dk} 



(26) 



Pr{g„GZ^fc} (27) 

JJf{9ml,gm2) dgmldgm2 



gmeOkHDi, 



pSel 
^k 
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Region Ok is the area within the border ek^Gek = 2 ^^*, which can be easily shown that, 
in terms of gi and 5^2, it is a hyperbola. 

C. Probability of rank failure 

As stated before, the destination after receiving M equations from the relays attempts to 
find and solve the two best independent equations, i.e. the two equations with the highest com- 
putation rates, to recover the both messages . Since the relays select their ECVs independently, 
the equations are not necessarily independent. If the destination fails to find two independent 
equations, a rank failure occurs and the destination cannot recover both transmitted messages. 
This event causes an outage as well. 

A rank failure occurs if all the M equations are linearly dependent. In a two-dimensional 
space, two vectors a and b are linearly dependent only when they are in the same direction, 
i.e. a = ab, a G M. From Lemma 3 and the fact that ECVs are integer vectors, two ECVs are 
linearly dependent only when they are equal. Therefore, a rank failure occurs at the destination 
if all the received equations are the same. Thus for CMF(ir), we can compute the probability 
of rank failure as 

K 

Pfaii = X] ^^ "L^m = efc , Vm = 1, . . . , M} 



fc=i 



^^f2^r{^.^ = e,r (28) 



k=l 
K 



E (pn 



k=l 

where P^^^ is defined in (20) and (a) follows from the fact that relays select their ECVs 
independently. 

It is useful to find a lower bound for the overall outage probability of the system. As stated 
above, a rank failure, with probability one causes an outage in the system. Hence, we can write 



psys _ p I n — P "i p^y^ 

^ out ^ fml ^ \^ ^ fail) --^ out\Nofail , , 

' (29) 

> Pfail 

where P^^l is the overall outage probability of the system and Plut\No fail ^^ ^^^ overall outage 
probability of the system conditioned on no rank failure. A direct result of the above equation 
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is that Pfaii lower bounds the overall outage probability of the system. 

It is noteworthy that (29) is true for the optimum CMF method as well (with a different 
probability of rank failure Pjaii)- Hence the probability of rank failure imposes a lower bound 
on the outage probability of the optimum CMF method as well. 

D. The outage probability of the system 

The destination receives M equations as well as their corresponding ECVs {s.m}m=i ^^'^ 
the computation rates {Rm\m=v Then, the receiver selects two independent equations with 
the highest rates. An outage occurs at the destination if either the minimum rate of these 
two selected equations is less than the given target rate Rt or the destination cannot find two 
independent ECVs among M received ECVs. Hence, the system outage probability can be 
expressed as 

PZ = Pr i . max . min {R, , R,) < rA . (30) 

I l^i,j^M, i^j, I 

\^ deti[ai,aj])jtO ) 

The probability in the above equation includes the event of rank failure in the destination, 
i.e. the case that destination cannot find two independent equations. In this case, the search 
space of maximization in (30) is an empty set and hence considered as an outage event. 

To compute the probability in (30) for the proposed scheme of CMF(K), we present a 
strategy for selecting two independent equations (we have assumed that L = 2) with the 
highest minimum rates, as follows: All the M received equations are divided into K sets, 
each having the same ECV. Please note that some of these sets can be empty. Then, the 
equation with the highest rate is selected from the non-empty sets. As a result, at most K 
independent equations, each from a non-empty set, are selected. Finally, the two equation with 
the highest rates among these independent equations are the desired ones. In the following, we 
calculate (30) based on this strategy. 

Define T^, A; = I, , K, as the set of received equations that their ECVs are e^. The size of 
the set Tfc is denoted by n^. From the definition of P^^' in (20), the probability of having the 
size vector of (tt-i, ^2, ■ ■ ■ , uk) is easily computed as 



Pr(ni,--- ,nK) 



'^ M ^ ^ 



ni,--- ,nK 



V 

13 



UiP^'T^ (31) 



k=l 



where 

I " \ 



Ml , , 



Thus, using the law of total probability, we drive 

Kll= E PrK'---'^^)-C*W..-,nK)' (33) 

(ni,---,nK) 
Vfc:Os£nfcs£M 

where P^^^\(n^ ■■■ n ) ^^ ^^^ outage probability given the size vector (ni, 77,2, ■ ■ ■ , riK). To compute 
the conditional outage , let K^ denotes the maximum computation rate of all equations in T^ 
. If Tfc is an empty set we set K^ = 0. Moreover, let the set {Rk,i, Rk,2,- ■■ ,Rk,nk} include 
all computation rates of equations in T^. Since all relays have the same conditions, all the 
computation rates are i.i.d., and we can calculate 

PT{R'k<Rt} 

= Pr {max {Rk,i, ■■■ , Rk,n^) < Rt} 

= Pi' s n ■^''''" < -^* r (34) 

= l[Pl{Rk,r<Rt} 

r=l 

= {Pout\ek} , 



where Pout\ek is given in (27). 

The conditional probability in (33) can be then computed as (35), where max(2) (■) repre- 
sents the second maximum operator and (a) follows from the fact that all the computation 
rates are i.i.d.. 
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(35) 



C'in,...,n.) = Pr {max(2) {R\, R'2, ■ ■ ■ , R'k) < Rt} 

= FTllf]R'k<RA\jl\jlR'k>Rt f] R'^<Rt 

[ \k=l J \k=l \ i=l,ji^k 

{K 1 ^ f ^ 

^R!k<Rt\ + Y.^AR!k>Rt n R!,<Rt 
k=l ) k=l I j=lj^k 

K K / K 

= l[PT{R'k<Rt} + J2hr{R'k>Rt}- n PHR'j<Rt} 

k=l k=l \ j=lJjLk 

= n iPout\.r + E f (1 - {Pou,.D ■ n iPou,.r 

k=l fc=l \ j=ljjtk 

Finally, by substitution of (31) and (35) in (33), the overall system outage probability is 
derived as (36), where P^^' and Pout\ek ^-^e given in (20) and (27), respectively. 

( M \ ^ 

ni,- ■■ ,nK / fc=i 



(ni,---,nK) \ 

n(^o-iej"+i:f(i-(^o«,ej"i- n iPou,.r)] (36) 

^fc=l fc=l \ 3=1 j^k J J 

V SIMULATION RESULTS 

In this section, computer simulations and analytical results are provided to study and 
compare optimum and simplified methods. In simulations, equal powers Pi = P2 = P and 
target rate Rt = 0.5 are assumed. 

Fig. 4 shows the probability of rank failure in the destination for optimum CMF, CMF(3), 
and CMF(5) methods versus the average SNR. Moreover, two cases of M = 2 and M = 6 
are compared. As it is observed, by increasing the number of relays (M), probability of rank 
failure decreases significantly. The reason is that by the increase of M, the destination receives 
more equations and hence it can find two independent equation with a higher probability. 
Generally, the optimum CMF shows a lower probability of rank failure compared with CMF(3) 
and CMF(5). Similarly, the probability of rank failure of CMF(5) is less than that of CMF(3). 
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This is due to the fact that in CMF(5), the search set has more elements than in CMF(3). 
Hence, the probabihty of selecting the same ECVs in different relays decreases. 

Fig. 5 compares the systems outage probability of optimum CMF, CMF(3), and CMF(5) 
versus average SNR in the cases of M = 2 and M = 6 relays, based on computer simulations. 
Analytical results are also provided for our proposed simplified schemes for the comparisons. 
First, as can be realized, simulation results for CMF(3) and CMF(5) well coincide with the 
related analytical results. From this figure, CMF(3) and CMF(5) outage curves provide upper 
bounds on optimum CMF outage curves, which are closely tight in SNRs less than 6 dB and 
16 dB, respectively. As it is observed, CMF(5) performs near the optimum CMF method up 
to a higher SNR threshold and provides a wider valid approximation region than CMF(3). 

Fig. 6 demonstrates the outage probability along with the probability of rank failure versus 
the average SNR, for the cases of optimum CMF and CMF(5). Here, M = 6 relays are assumed. 
As mentioned in Subsection IV-C and can be observed from this figure, probability of rank 
failure lower bounds the outage probability in the both cases. The lower bound is closely 
tight in SNRs more than 12 dB. Therefore, rank failure is a major bottleneck in the system 
and degrades the performance considerably, even in the optimum CMF method. From these 
results, by using a higher number of relays (M > L), the rank failure probability reduces which 
leads to a lower outage probability. 

Fig. 7 illustrates the effect of CEE on the optimum and simplified CMF methods. In this 
figure, the outage probabilities of optimum CMF and CMF(5) versus the average SNR are 
shown for different values of CEE variance a^ through computer simulations. The number of 
relays (M) is assumed to be 6. In the simulations, the true channel gain hmi and the estimated 
channel gain hmi are computed as 

hml = \lml\ , (37) 

hml = \lml + <yeeml\ , (38) 

where 7^; and Cmi, for / = l,2,m = l,---, M, are i.i.d circularly-symmetric complex normal 
random variables with unit variance. CEE results in a noisy matrix Gm in (6) and (17) and 
hence, a suboptimum ECV may be selected as the solution of these optimizations. Selection of a 
suboptimum ECV causes a rate loss. It can be found from this figure that the optimum CMF is 
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more sensitive to CEE than the simphfied CMF, and even in some cases shows shghtly inferior 
performance than the simphfied CMF. The reason is that optimum CMF has a larger search 
space, especiahy in high SNR values, than the simplified CMF, which has a small and fixed 
search space. Hence, when G^ is noisy due to CEE, it is more probable that a suboptimum 
ECV is selected as the solution of (6) than (17). 

VI CONCLUSION 

In this paper, we have proposed a simplified CMF method, called CMF{K). Through ana- 
lytical and simulation evaluations, it has been demonstrated that the CMF(i^) presents nearly 
the same performance as optimum CMF in low SNR regimes, where the SNR is below a cer- 
tain threshold. Higher values of K increase this threshold and extend the valid approximation 
region. The exact outage probability of CMF(K) is derived. Here, we have considered two 
sources and real channels. However, the same framework can be applied to extend the results 
to a greater number of sources or complex cases. Our results indicates that the rank failure of 
the received equations at the destination is a dominate performance degradation in term of the 
overall system outage probability in both optimum and our proposed simplified CMF, which 
its effect can be reduced by having higher number of relays compared to the number of sources. 
Finally, computer simulations showed that simplified CMF, due to its small and fixed search 
space, is more robust against channel estimation error than optimum CMF for the examples 
considered. 
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Fig. 2: Probability of selecting different ECVs in a relay for optimum CMF method versus 
average SNR (Pi = Pg = P)- 
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Fig. 4: Probability of rank failure in destination for optimum CMF, CMF(3), and CMF(5) 
versus average SNR (M = 2, 6). 
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Fig. 5: System outage probability for optimum CMF, CMF(3), and CMF(5) versus average 
SNR, (M = 2,6). 
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Fig. 6: Probabilities of outage and rank failure for optimum CMF and CMF(5) versus average 
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